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Abstract 

We prove almost sure Euler hydrodynamics for a large class of attrac- 
tive particle systems on Z starting from an arbitrary initial profile. 
We generalize earlier works by Seppalainen (1999) and Andjel et al. 
(2004). Our constructive approach requires new ideas since the subad- 
ditive ergodic theorem (central to previous works) is no longer effective 
in our setting. 



1 Introduction 

Hydrodynamic limit ([131 EH [211 HQ]) is a law of large numbers for the time 
evolution (usually described by a limiting PDE, called the hydrodynamic 
equation) of empirical density fields in interacting particle systems {IPS). In 
most known results, only a weak law of large numbers is established. In this 
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description one need not have an explicit construction of the dynamics: the 
hmit is shown in probabihty with respect to the law of the process, which is 
characterized in an abstract way by its Markov generator and Hille-Yosida's 
theorem ([29]). Nevertheless, when simulating particle systems, one naturally 
uses a pathwise construction of the process on a Poisson space-time random 
graph (the so-called "graphical construction"). In this description the dy- 
namics is deterministically driven by a random space-time measure which 
tells when and where the configuration has a chance of being modified. It 
is of special interest to show that the hydrodynamic limit holds for almost 
every realization of the space-time measure, as this means a single simulation 
is enough to approximate solutions of the limiting PDE. In a sense this is 
comparable to the interest of proving strong (rather than weak) consistency 
for statistical estimators, by which one knows that a single path observation 
is enough to estimate parameters. 

Most usual IPS can be divided into two groups, diffusive and hyperbolic. 
In the first group, which contains for instance the symmetric or mean-zero 
asymmetric exclusion process, the macroscopic^microscopic space-time scal- 
ing is {x,t) I— i> {Nx,N'^t) with N oo, and the limiting PDE is a diffusive 
equation. In the second group, which contains for instance the nonzero mean 
asymmetric exclusion process, the scaling is (x, t) ^ {Nx, Nt), and the limit- 
ing PDE is of Euler type. In both groups this PDE often exhibits nonlinearity, 
either via the diffusion coefficient in the first group, or via the flux function 
in the second one. This raises special difficulties in the hyperbolic case, due 
to shocks and non-uniqueness for the solution of the PDE, in which case the 
natural problem is to establish convergence to the so-called "entropy solu- 
tion" ([SZl)- In the diffusive class it is not so necessary to specifically estabhsh 
strong laws of large numbers, because one has a fairly robust method ([25]) 
to establish large deviations from the hydrodynamic limit. Large deviation 
upper bound and Borel Cantelli's lemma imply a strong law of large numbers 
as long as the large deviation functional is lower-semicontinuous and has a 
single zero. The situation is quite different in the hyperbolic class, where 
large deviation principles are much more difficult to obtain. So far only the 
remarkable result of Jensen and Varadhan ([21] and [H]) is available, and 
it applies only to the one-dimensional totally asymmetric simple exclusion 
process. Besides, the fact that the resulting large deviation functional has 
a single zero is not at all obvious: it follows only from recent and refined 
work on conservation laws ([12j). An even more difficult situation occurs for 
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particle systems with nonconvex and nonconcave flux function, for which the 
Jensen- Varadhan large deviation functional does not have a unique zero, due 
to the existence of non-entropic solutions satisfying a single entropy inequal- 
ity (|20]). In this case a more complicated rate functional can actually be 
conjectured from |H [3T]. 

The derivation of hyperbolic equations as hydrodynamic limits began with 
the seminal paper [3^, which established a strong law of large numbers for 
the totally asymmetric simple exclusion process on Z, starting with I's to 
the left of the origin and O's to the right. This result was extended by [5] 
and [T] to nonzero mean exclusion process starting from product Bernoulli 
distributions with arbitrary densities A to the left and p to the right (the so- 
called "Riemann" initial condition). The Bernoulli distribution at time is 
related to the fact that uniform Bernoulli measures are invariant for the pro- 
cess. For the one-dimensional totally asymmetric iT-exclusion process, which 
does not have explicit invariant measures, a strong hydrodjTiamic limit was 
established in starting from arbitrary initial profiles, by means of the 
so-called "variational coupling". These are the only strong laws available 
so far. A common feature of these works is the use of subadditive ergodic 
theorem to exhibit some a.s. limit, which is then identified by additional 
arguments. 

On the other hand, many weak laws of large numbers were established for 
attractive particle systems. A first series of results treated systems with 
product invariant measures and product initial distributions. In [2], for a 
particular zero-range model, a weak law was deduced from conservation of 
local equilibrium under Riemann initial condition. It was then extended in 
[3] to the misanthrope's process of [TT] under an additional convexity as- 
sumption on the flux function. These were substantially generalized (using 
Kruzkov's entropy inequalities, see [2^) in [32] to multidimensional attrac- 
tive systems with arbitrary Cauchy data, without any convexity requirement 
on the flux. For systems with unknown invariant measures, the result of 
[35i| was later extended to other models, though only through weak laws. In 
[33] , using semigroup point of view, hydrodynamic limit was established for 
the one-dimensional nearest-neighbor i^-exclusion process. In [7] we studied 
fairly general attractive systems, employing a constructive approach we had 
initiated in |6]. This method is based on an approximation scheme (to go 
from Riemann to Cauchy initial data) and control of some distance between 
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the particle system and the entropy solution to the hydrodynamic equation. 

In this paper we prove a strong law of large numbers, starting from arbi- 
trary initial profiles, for finite-range attractive particle systems on Z with 
bounded occupation number. We need no assumption on the flux, invari- 
ant measures or microscopic structure of initial distributions. This includes 
finite-range exclusion and ii'-exclusion processes, more general misanthrope- 
type models and fc-step exclusion processes. We proceed in the constructive 
spirit of [H [7], which is adapted to a.s. convergence, though this requires a 
novel approach for the Riemann part, and new error analysis for the Cauchy 
part. The implementation of an approximation scheme involves Riemann 
hydrodynamics from any space-time shifted initial position, which cannot be 
obtained through subadditive ergodic theorem, and thus prevents us from 
using the a.s. result of [Ij. This can be explained as follows (see Appendix 
|B]for more details). Assume that on a probability space we have a shift op- 
erator 9 that preserves the probability measure, and a sequence of 
real-valued random variables converging a.s. to a constant x. Then we can 
conclude that the sequence '■= ° d"')n converges in probability to x, 
but not necessarily almost surely. Indeed, for fixed n, Yn has the same distri- 
bution as Xn, but the sequence (Yn)n need not have the same distribution as 
(X„)„. So the derivation of a.s. convergence for Yn is a case by case problem 
depending on how it was obtained for X„. In particular, if convergence of 
Xn was established from the subadditivity property 

Xn+m < Xn + X^ O 0^ 

this property is no longer satisfied by Yn- In contrast, if convergence of Xn 
was established by large deviation estimates for X„, these estimates carry 
over to Yn, thus also implying a.s. convergence of {Yn)n- In our context the 
random variable X„ is a current in the Riemann setting. In order to han- 
dle the Cauchy problem we have to establish a.s. convergence for a shifted 
version Yn of this current. We derive asymptotics for Yn by means of large 
deviation arguments. 

The paper is organized as follows. In Section [2], we define the model, give its 
graphical representation, its monotonicity properties, and state strong hy- 
drodynamics. In Section [31 we derive almost sure Riemann hydrodynamics. 
In Section m we prove the result starting from any initial profile. 



4 



2 Notation and results 



Throughout this paper N = {1, 2, ...} will denote the set of natural numbers, 
Z+ = {0, 1, 2, ...} the set of non-negative integers, and R+* = R+ - {0} the 
set of positive real numbers. We consider particle systems on Z with at 
most K particles per site, e N. Thus the state space of the process is 
X = {0, 1, ■ ■ ■ , K}"^, which we endow with the product topology, that makes 
X a compact metrisable space. A function defined on X is called local if it 
depends on the variable r] & X. only through {ri{x),x e A) for some finite 
subset A of Z. We denote by t^, either the spatial translation operator on the 
real line for a; G M, defined by r^y = x -\- y, or its restriction to Z for a; £ Z. 
By extension, we set Txf = /or^ if / is a function defined on M; r^rj = rjoTx, 
for a; e Z, if 7] G X is a particle configuration (that is a particular function 
on Z); TxH = no r^^ if /x is a measure on R or on X. We let A^+(]R) denote 
the set of positive measures on R equipped with the metrizable topology 
of vague convergence, defined by convergence on continuous test functions 
with compact support. The set of probability measures on X is denoted by 
P(X). If T] is an X- valued random variable and u e 7^(X), we write ~ 
to specify that rj has distribution u. The notation where / is a real- 

valued function and u e ■P(X), will be an alternative for J^fdu. We say 
a sequence G N) of probability measures on X converges weakly to 

some z/ G V(X.), if and only if as n — oo for every continuous 

function / on X. The topology of weak convergence is metrizable and makes 
V(X.) compact. 

2.1 The system and its graphical construction 

We consider Feller processes on X whose transitions consist of particle jumps 
encoded by the Markov generator 

Lfiv) -J2p(y- ^)Kv{x),v{y)) [f iv^'") - fiv)] (1) 

x,yeZ 

for any local function /, where 7]^'^ denotes the new state after a particle has 
jumped from x to y (that is ri^''^{x) = rj^x) — 1, ri^'y{y) = rjiy) + 1, rj^'y^z) = 
r){z) otherwise), p is the particles' jump kernel, that is a probability distri- 
bution on Z, and b : Z+ x Z+ — > 1R+ is the jump rate. We assume that p 
and b satisfy : 
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(Al) The semigroup of Z generated by the support of p is Z itself 
(irreducibihty) ; 

(A2) p has a finite first moment, that is J2zez \^\pi^) < +00; 
(A3) 6(0, .) = 0, 6(., K) = (no more than K particles per site), and 
6(1, is: - 1) > 0; 

(A4) 6 is nondecreasing (nonincreasing) in its first (second) argument. 

We denote by {S{t),t G M^) the semigroup generated by L. Without addi- 
tional algebraic relations satisfied by 6 (see [H]), the system in general has 
no explicit invariant measures. This is the case even for K- exclusion process 
with K >2 (see [3^), where b{n,m) = l{n>o,m<K}- 

Remark 2.1 For the sake of simplicity, we restrict our attention to processes 
of "misanthrope type " , for which particles leave more likely crowded sites for 
less occupied ones (cf. ITT^)- But our method can he applied to a much 
larger class of attractive models that can he constructed through a graphical 
representation, such as k-step exclusion (see fTl^ and J^). 



We now describe the graphical construction of the system given by ([T]), which 
uses the Harris representation (^I9j, ^ p. 172], [8, p. 119], [301 P- 215]); 
see for instance [T] for details and justifications. This enables us to define 
the evolution from different initial configurations simultaneously on the same 
probability space. We consider the probability space (f2, JF, P) of measures 
on M+ X Z2 X [0, 1] of the form 



uj{dt, dx, dz, du) = ^ 5{t„,x„ 



,Um.) 



mGN 



where denotes Dirac measure, and {tm,Xm, Zm,Um)m>o are pairwise dis- 
tinct and form a locally finite set. The u-field J-" is generated by the mappings 
u ^— s> uj{S) for Borel sets S. The probability measure P on i7 is the one that 
makes u a Poisson process with intensity 

m{dt,dx,dz,du) = ||6||ooAM+((it) x Xi,{dx) x p{dz) x X[Q ij(du) 

where A denotes either the Lebesgue or the counting measure. We denote by 
IE the corresponding expectation. Thanks to assumption (A2), for P-a.e. u, 
there exists a unique mapping 

(r^o, t) G X X M+ ^ = 77^(7/0, cu) G X (2) 
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satisfying: (a) t i-^ rjt{rjQ,u) is right-continuous; (b) rjQ{rjQ,uj) = rjo; (c) for 
teR+, ix,z) EZ^,rit = ri^r'+' if 

3^6 [0,1]: u;{(t,x,z,M)} = landM<^^^^^^^^^(f^^^ (3) 



and (d) for all s, t G M+* and a; G Z, 

uj{[s,t] X Z^x (0, 1)} = ^ Vi; G [s,t],r7„(x) = r7,(x) (4) 

where 

:= {{y, z) e Z"^ : y = X or y + z = x] 

In short, ([3]) tells how the state of the system can be modified by an "a;- 
event", and (jl]) says that the system cannot be modified outside a;-events. 
This defines a Feller process with generator ([1]): that is for any t G M"*" and 
/ G C(X) (the set of continuous functions on X), S(t)f G C(X) where 

sit)fivo) = nfivtim,uj))]. 

An equivalent formulation is the following. For each {x, z) G Z^, let {T^'^, n > 
1} be the arrival times of mutually independent rate ||6||ooP(-2) Poisson pro- 
cesses, let {U^'^,n > 1} be mutually independent (and independent of the 
Poisson processes) random variables, uniform on [0, 1]. At time t = T^'^, the 

b{r]t-{x),r]t-{x + z)) 

It- " " " 

unchanged otherwise. 



configuration r]t- becomes rj^f ^ if f/^''^ < ^ " ^ ^ — ' — —, and stays 



Remark 2.2 For the K-exclusion process, b takes its values in {0, 1}, the 
probability space can be reduced to measures u{dt,dx,dz) on IR+ x Z^, and 
^ to r]t = ri^f^^ if ?7i- (x) > and r]t- {x + z) < K. One recovers exactly 
the graphical construction presented in e.g. fW^ or fJ^. 

One may further introduce an "initial" probability space (f2o, .Fq, Pq), large 
enough to construct random initial configurations r/o = Voi^o) for cuq G Qq. 
The general process with random initial configurations is constructed on the 
enlarged space {Q, = ^Iq x il, T = cr(J?^o x JF), P = Pq ® P) by setting 

r]t{uj) = r]t{r]o{uJo),uj) 

for uj = (ljq, uj) G Q. If TjQ has distribution /iq, then the process thus con- 
structed is Feller with generator ([T]) and initial distribution fiQ. By a coupling 
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of two systems, we mean a process {rjt, Ct)t>o defined on Q, wliere eacli com- 
ponent evolves according to and tlie random variables rjo and Co are 
defined simultaneously on Qq. 

We define on fl the space-time shift 0^0^: for any u E Q, for any {t, x, z, u) 

{t, X, z, u) G OxQ^to^ if and only if {to + t,XQ + x, z,u) E iv 

where (t,x,z,u) G u means u{{t,x, z,u)} = 1. By its very definition, the 
mapping introduced in ([2]) enjoys the following properties, for all s,t > 0, 
X e Z and {rj, G X x fi: 

VsiVtiv, %,t^) = Vt+siv, ^) (5) 
which implies Markov property, and 

^xVtiV, ^) = VtiTxV, dxfi^) (6) 

which implies that S{t) and r^, commute. 



2.2 Main result 

We give here a precise statement of strong hydrodynamics. Let G N be 
the scaling parameter for the hydrodynamic limit, that is the inverse of the 
macroscopic distance between two consecutive sites. The empirical measure 
of a configuration rj viewed on scale is given by 

a^{T^){dx) = N-^Y^T^{y)5y/N{dx) G M+{M) 
We now state our main result. 

Theorem 2.1 Assume p{.) is finite range, that is there exists M > such 
that p{x) = for all \x\ > M. Let {rj^ , N E N) be a sequence of X- 
valued random variables on VLq. Assume there exists a measurable [0,/^]- 
valued profile Uq{.) on M such that 

lim {riQ){dx) = Uo{.)dx, Pq-o-s. (7) 

TV— +00 

that is, 

lim / ijj{x)a'^{riQ){dx) = / ^lj{x)uo{x)dx, Pq-o-s. (8) 
Jr J 
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for every continuous function ip on M. with compact support. Let {x, t) 
u{x,t) denote the unique entropy solution to the scalar conservation law 

dtU + d^[G{u)]=0 (9) 

with initial condition uq, where G is a Lip schitz- continuous flux function (de- 
fined in ( f7^) below) determined by p{.) and b{.,.). Then, with JP -probability 
one, the convergence 

lim a^{r]m{Vo {^o),^^)){dx) = u{.,t)dx (10) 

holds uniformly on all bounded time intervals. That is, for every continuous 
function iponW with compact support, the convergence 

lim (11) 

holds uniformly on all bounded time intervals. 

As a byproduct we derive in Appendix |X] the main result of [7j: 

Corollary 2.1 The strong law of large numbers fT^) implies the weak law 
of large numbers established in [7f. 

We recall from [71 pp. 1346-1347 and Lemma 4.1] the definition of the Lipschitz- 
continuous macroscopic flux function G. For p G 7^ C [0,K] (the closed set 
of allowed densities for the system, defined in Proposition 12.11 below) let 



G{p) = 



J^zp(z)6(r/(0),r/(z)) 



:i2) 



this represents the expectation, under the shift invariant equilibrium measure 
with density p (see Proposition 12.11) . of the microscopic current through site 
0. On the complement of 7^, which is at most a countable union of disjoint 
open intervals, G is interpolated linearly. A Lipschitz constant of G is 
determined by the rates 6, p in ([1]) : 

V = 2B'^\z\p{z), with 

B = sup {b{a,k)—b{a,k-\-l),b{k-\-l,a)—b{k,a)} 

0<a<K,0<k<K 
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This constant will be a key tool for the finite propagation property of the 
particle system (see Proposition 14. II below). 

In [71 Section 2.2] we discussed the different definitions of entropy solutions 
(loosely speaking : the physically relevant solutions) to an equation such as 
(jni), and the ways to prove their existence and uniqueness. Therefore we just 
briefly recall the deflnition of entropy solutions based on shock admissibility 
conditions (Olemik's entropy condition), valid only for solutions with locally 
bounded space-time variation ([12]). A weak solution to (Q is an entropy 
solution if and only if, for a.e. t > 0, all discontinuities of u{.,t) are entropy 
shocks. A discontinuity («",«"•"), with := u{x ± 0, t), is called an entropy 
shock, if and only if: 

The chord of the graph of G between u~ and lies 
below the graph if < u^, above the graph if < . 

In the above condition, "below" or "above" is meant in wide sense, that is 
the graph and chord may coincide at some points between u~ and 
If the initial datum Uq{.) has locally bounded space- variation, there exists a 
unique entropy solution to ([2]), within the class of functions of locally bounded 
space-time variation ([12]). 

2.3 Monotonicity and invariant measures 

Assumption (A^) implies the following monotonicity property, crucial in our 
approach. For the coordinatewise partial order on X, deflned by ?7 < C if and 
only if ri{x) < (^{x) for every a; G Z, we have 

(?7o,t) ^ rit{riQ,uj) is nondecreasing w.r.t. rjo (13) 

for every uj such that this mapping is well deflned. The partial order on X 
induces a partial stochastic order on 'P(X); namely, for /ii,/i2 G '^(X), we 
write /ii < fi2 if the following equivalent conditions hold (see e.g. [29], [39]): 

i) For every non- decreasing nonnegative function / on X, /ii(/) < /^2(/)- 

ii) There exists a coupling measure /I on X = X x X with marginals ni and 
fi2, such that /i{(?7, ^) : r] < = 1. 

It follows from this and f[T^ that 

/^i < /i2 ^ Vt G M+, fiiS{t) < fi2S{t) (14) 
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Either property (fT^ or (HM is usually called attractiveness. 



Let X and S denote respectively the set of invariant probability measures 
for L, and the set of shift-invariant probability measures on X. We derived 
in [71 Proposition 3.1], that 

Proposition 2.1 

(J^5), = {z/^pG7^} (15) 

with TZ a closed subset of [0, K] containing and K , and a shift-invariant 
measure such that i^P[ri{0)] = p. (The index e denotes extremal elements.) 
The measures are stochastically ordered: 

p<p'^uP< yP' (16) 

Moreover we have, quoting [331 Lemma 4.5]: 

Proposition 2.2 The measure vP has a.s. density p, that is 

1 ' 

lim — > nix) = p, fP — a.s. 

x=—l 

By [22^ Theorem 6], f|T6l) implies existence of a probability space on which 
one can define random variables satisfying 

r]P ~ lyP (17) 

and, with probability one, 

r]P < r]p', Vp, p' en such that p < p' (18) 

Proceeding as in the proof of [23, Theorem 7], one can also require (but 
we shall not use this property) the joint distribution of {rjP : p G 7?.) to be 
invariant by the spatial shift t^. In the special case where uP are product 
measures, that is when the function &(.,.) satisfies assumptions of [TT], such 
a family of random variables can be constructed explicitely: if {Ux)^^^^ is a 
family of i.i.d. random variables uniformly distributed on (0, 1), one defines 

V^'ix) = F-\U.,) (19) 
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where Fp is the cumulative distribution function of the single site marginal 
distribution of u^. We will assume without loss of generality (by proper en- 
largement) that the "initial" probability space is large enough to define a 
family of random variables satisfying f|T7|) - f|T8|) . 

An important result for our approach is a space-time ergodic theorem for 
particle systems mentioned in [Mj, which we state here in a general form, 
and prove in Appendix O 

Proposition 2.3 Let {rit)t>o be a Feller process on X with a translation 
invariant generator L, that is 

TiLr.i = L (20) 

Assume further that 

/i G (Xl n S)e 

where Xl denotes the set of invariant measures for L. Then, for any local 
function f on X, and any a > 

fc^X g^-^''")* = y^* = £Ssj/„ (21) 

a.s. with respect to the law of the process with initial distribution ^. 
Remark 2.3 It follows from [21\} that, more generally, 

1 f^^ f 

}'SL (i-a)(d-c)p L ^ = (22) 

^ i(zzn[ci,di] 

for every < a < b and c < d in M., as can be seen by decomposing the 
space-time rectangle [al, hi] x [c/, dl] into rectangles containing the origin. 



3 Almost sure Riemann hydrodynamics 

By definition, the Riemann problem with values A, p G [0, K] for (j9]) is the 
Cauchy problem for the particular initial condition 

Rlpix) = Xl{,^o} + pM.m (23) 
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The entropy solution for this Cauchy datum will be denoted in the sequel by 
Rx^p{x,t). In this section, we derive the corresponding almost sure hydrody- 
namic limit when \, p & TZ. We will use the following variational represen- 
tation of the Riemann problem. We henceforth assume X < p (for the case 
X > p, replace minimum with maximum below and make subsequent changes 
everywhere in the section). 

Proposition 3.1 (I?, Proposition 4-1])- Let A,p G [Oj-ft'], X < p. 
i) There is a countable set S/o^ C [A,p] (depending only on the differentia- 
bility properties of the convex envelope of G on [X, p]) such that, for every 
i; G M \ T^iow, G{-) — V- achieves its minimum over [A,p] at a unique point 
hciy). Then R\^p{x,t) = hc{x/t) whenever x/t ^ Sjow 

a) Suppose X,p Then the previous minimum is unchanged if restricted 

to [X, p] mZ. As a result, the Riemann entropy solution is a.e. TZ-valued. 

Remark 3.1 Property ii) holds if we replace TZ by any closed subset of [0, K] 
on the complement of which G is affine. We stated it for TZ because it is the 
set of densities relevant to the particle system. 

To state Riemann hydrodynamics, we define a particular initial distribution 
for the particle system. We introduce a transformation T : ^ X by 

We define z/^'^ as the distribution of T{ri^, rjp) =: rj^'f, and z/^'^ as the coupling 
distribution of {r]'^,r]P). Note that, by ffTS]) . 

z>^'^{(r7,OeX2: r/<^} = l (24) 

The measure v^'^ is non-explicit unless we are in the special case of [TT] where 
the i/P are product, one can use ffT^ . and z/^'^ is itself a product measure. In 
all cases, v'^'P enjoys the properties: 

PI) Negative (nonnegative) sites are distributed as under (z/''); 

P2) Tiu^^P > v^'P {nv^^P < u^'P) if A < p (A > p); 

P3) u'^'P is stochastically increasing with respect to A and p. 

Let us also define an extended shift 6' on the compound probability space 
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Q' = X Q. This is a particular case of Q when the "initial" probability 
space Qq is the set of coupled particle configurations (rj,^). Let 

uj' = {v,^,uj) (25) 
denote a generic element of this space. We set 

= i'TxVtiv^ r^Vti^, O^^tuj) (26) 
We can now state and prove the main results of this section. 
Proposition 3.2 Set 

Afr iuj'):= J2 Vmv,0,u^)ix) (27) 

[vt]<x<lwt] 

Then, for every t > 0, a G M^, /? G M and f , w G M \ 
hm V''"o^U],„.(^') 

= [GiRx,piv, 1)) - V Rx,p{v, 1)] - [G(i?A,p(w;, 1)) - w Rx,p{w, 1)] (28) 
v^'p (g) P-a.s. 

Remark 3.2 T/izs result is an almost sure version of [3, Lemma 3.2], where 
the limit of the corresponding expectation was derived. 

Corollary 3.1 Set 

f3f{u;'){dx) := a'^iMTiv^O^^mdx) (29) 
(i) For every t > 0, Sq > and xq G M, we have the z/'^'^® P-a.s. convergence 
lim PNti^[Nxo].Nso^')idx) =Rx,p{.,t)dx 

TV — ^oo 

(a) In particular, for an initial sequence {ri^)N such that rj^ = rj^'P for 
every G N, the conclusion of Theorem \2.1\ holds without the finite-range 
assumption onp{.). 
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For the asymmetric exclusion process, fl] proved a statement equivalent to 
the particular case a = /3 = of Proposition I3.2[ Their argument, which is 
a correction of [5], is based on subadditivity. As will appear in Section HI we 
do need to consider nonzero a and f3 in order to prove a.s. hydrodynamics 
for general (non-Riemann) Cauchy datum. Using arguments similar to those 
in [T] , it would be possible to prove Proposition 13.21 with a = (3 = for our 
model ([T]). However this no longer works for nonzero a and (3 (see Appendix 
[B|) . Therefore we construct a new approach for a.s. Riemann hydrodynam- 
ics, that does not use subadditivity. 

To prove Proposition 13. 2[ we first rewrite (in Subsection 13.11) the quantity 
JVf'^ in terms of particle currents for which we then state (in Subsection 13.21) 
a series of lemmas (proven in Subsection 13.41) . and finally obtain the desired 
limits for the currents, by deriving upper and lower bounds. For one bound 
(the upper bound if A < p or lower bound if A > p), we derive an "almost-sure 
proof" inspired by the ideas of |3j and their extension in [7]. For the other 
bound we use new ideas based on large deviations of the empirical measure. 

3.1 Currents 

Let us define particle currents in a system {rjt)t>Q governed by ([3])-(jl]). Let 
X. = (xt, t > 0) be a Z- valued cadlag path with \xt — Xt- 1 < 1. In the sequel 
this path will be either deterministic, or a random path independent of the 
Poisson measure u. We define the particle current as seen by an observer 
travelling along this path. We first consider a semi-infinite system, that is 
with J2x>o^o(-^) ^ +oo: in this case, we set 



where rjt^riQ, u) is the mapping introduced in ([2]). In the sequel we shall most 
often omit u and write rit for 774(770,0;) when this creates no ambiguity. We 
have 



(30) 



y>xt y>xo 



(31) 
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where 



V^t" iVo,^) = uj <^{s,y,z,u) : < s <t, y < Xs < y + z, 

. Kvs-{y),Vs-{y + z)) 

Vt"~{m,i^) = l^{s,y,z,u) : < s <t, y + z < Xs < y, 

. Kvs-{y),ris-{y + z)) 
u < 



(32) 



count the number of rightward/leftward crossings due to particle jumps, and 



ip1-{r]Q,uj) = - I ris-{xs\/ Xs-)dxs (33) 
J[o,t] 

is the current due to the self-motion of the observer. For an infinite system, 
we may still define ip^-{r]o,uj) by equations (!3T!) to (!33|) . We shall use the 
notation Lfl in the particular case Xt = [vt]. The following identities are 
immediate from (1301) in the semi-infinite case, and extend to the infinite 
case: 

IftiVo,^) - ffiVo,^)] < K {\xt - yt\ + \xo - yo\) (34) 

[wt] 

Vt{m,^){x) = <^t{vo,i^) - <^T{Vo,^) (35) 

x=l+[vt] 

Following (|35|1 . the quantity Mt'^iuj') defined in fl27|) can be written as 

ur{u;') = <Pi{^')-<pr{^') (36) 

where we define the current 

0^(a.'):=^l'(T(r/,O,^) 

Notice that 

The proof of the existence of the limit in Proposition 13.21 is thus reduced to 
lim (d'm at^') exists u^'^ ® P-a.s. (37) 
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3.2 Lemmas 

We fix a G and /9 G M. The first lemma deals with equilibrium processes. 

Lemma 3.1 For all r E [\, p] n TZ, <^ e X, v E R \ T^iou,, 

lim t- o OL. uj) = G{r) - vr, z7^'^ ® P-a.s. 

The second lemma relates the current of the process under study with equi- 
librium currents; here it plays the role of Lemma 3.3 in [3j. It is connected 
to the "finite propagation property" of the particle model (see Lemma : 

Lemma 3.2 There exist v and v (depending on b and p) such that we have, 
iy^'P (g) P-a.s., 

\\m [t- ° uj) - t-m o Q[^,^^^,{i, ^)] = 0, (38) 

jor all V > V, and 

\\m [t- ° o'mUv. e, uj) - 1- vr ° V, ^)] = 0, (39) 

for all V < V. 

For the next lemmas we need some more notation and definitions. Let f G M. 
We consider a probability space , whose generic element is denoted by 
on which is defined a Poisson process A^^^ = Nt{u~^) with intensity |f |. We 
denote by P"*" the associated probability. We set 

xl{u+) := {sgn{v)) [N^t+,{uj+) - N^t{uj+)] (40) 
Tf,{r]o,uj,uj'^) := T^i{u^+)ris{riQ,uj) (41) 

Thus {jfs)s>o is a Feller process with generator 

L, = L + 5,, 5,/(C) = [/(rsgn(.)C) - /(C)] (42) 

(for / local and C G X) for which the set of local functions is a core, as it is 
known to be ([29]) for L. We denote by the set of invariant measures for 
Ly. Since any translation invariant measure on X is stationary for the pure 
shift generator 5"^, we have 

ir\S = z,r\S (43) 
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It can be shown (see [Ul Theorem 3.1] and [121 Corollary 9.6]) that d S 
for f 7^ 0, which implies 

but we shall not use this fact here. Define the time empirical measure 

mi(cj',U;+) := t^W %t(T(T,,0,a;,a;+)C^S (44) 

Jo 

and space-time empirical measure (where £ > 0) by 

mt,e{oj\uj+):=\Zr\[-et,et]\-^ ^ r,mi(cu', o;+) (45) 

x&Z: \x\<et 

We introduce the set 

Mx,p ■■= G P(X) : z/^ < /i < z/^} 
Notice that this is a closed (thus compact) subset of the compact space 'P(X). 

Lemma 3.3 (i) With u^^'^^JP^JP^ -probability one, every subsequential limit 
(^f ^t,e{(^[i3t] at^' '^^^ as t ^ oo Ues in n S n Aix^p = X nS H A^A,p- 
(a) XnSnM.x,p is the set of probability measures v of the form = f i'^'j{dr), 
where is a probability measure supported on Tin [A,p]. 

The proof of Lemma 13.31 will be based on the following large deviation result 
in the spirit of [T^. 

Lemma 3.4 (i) The functional defined on V(X.) by 

V^{lj) := sup - / -^{ri)dfi{r]) (46) 

/local J eJ 

is nonnegative, lower-semicontinuous, andV~^{0) =Xv. 

(a) Let ^ be a Markov process with generator and distribution denoted by 
P. Define the empirical measures 

Ml):=t-' f 6^ds, 7Tt,e:=\y^n[-et,et]\-' ^ r^vr^ (47) 

° x&n[-et,et] 

where e > 0. Then, for every closed subset F ofV(X.) and every t > 0, 

limsup t-^ logP (rrt^eii) e f) < - inf (48) 
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3.3 Proofs of Lemma 13.11 Proposition 13.21 and Corol- 
lary D 



Proof of Proposition U7B . We denote by uj' = {ri,C,,uj) a generic element of 
Q'. We will establish the following limits: first, 

M [G{r)-vr] < liminf rV* o 

< limSUp 0^'/3t],at(^') 

< sup [G{r)-vr], z/^'^ ® P-a.s. (49) 

r€[x,p]nn 



and then 



limsup t-Vr o O'wt] ati^') < inf [G{r) - vr], u^'P ® P-a.s. (50) 

which will imply the result, when combined with Proposition 13. II and the ex- 
pression (136!) of Mt'^iuj'). Though only the first inequality in fH9l) (together 
with (!50|) ) seems relevant to Proposition 13. 2[ we will need the whole set of 
inequalities: Indeed, writing (H^ for A = p = r G 7^ proves the equilibrium 
result of Lemma 13.11 

To obtain the bounds in fj49l) we proceed as follows. First we replace the 
deterministic path vt in the current (pl by x\{(jJ^). Then we consider a spa- 
tial average of v'^ *''^ '^^^ for x G [—et.et]^ and we introduce, for C G X, the 
martingale M^''"{C,,uj^uj^) associated to v^^ '"'^ ''^^{C,uj) (see below (ITOI) ). An 
exponential bound on the martingale reduces the derivation of fH9l) to bounds 
(deduced thanks to Lemma [373|l on 

(see below (Mli)). where [f{ri) — 11^7(1)] corresponds to the compensator of 
(a; )+^(^^^|^^ '''(^, c<j, cu"''). The bound (15U]) relies on Lemmas 13.11 and 

13.21 combined with the monotonicity of the process. 



Step one: proof of (J^- We have 



t-'<Pt ° ^Uati^') = t-Mi^at, 0mM^) (51) 
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where the configuration 

■^at = ^ativ, ^, ^) ■= T (r[/3i]?7„t(r7, u), T^pt]riat{i, ^)) (52) 

depends only on the restriction of uj to [0, at] x Z. Thus, since is a Poisson 
measure under P, 0[jit],at^ is independent of Watijl^ ^) ^) under u^'P ® P, and 

under v^'P ® P C?) P^, vOat is independent of {6]^j3t]^oit^,uj^) (53) 

Define 

y£l:\y\<et 

for every {(,uj,uj^) G X x ri x ri+, with a;*(cj+) given by By (IMI) . 

|^;'(C,cu) - V^r(C,^,^+)| < K {2et+\x\{uj^)-vt\) 

Since t~^a;J(ci;+) — v with P^-probabihty one, the proof of (H9|) can be 
reduced to that of the same inequahties with the l.h.s. of (15T]) replaced by 

t"Vr^(Wai,6'[/3t],atU;,CJ+) (54) 

and v^^P <S) P replaced by u^'P (g) P (g) P+. Let f{r]) := f^{r]) - f'{T]), with 
f~^iv) = Y P{z)b{r]{y),r]{y + z)) 

y^zaT,: y+z<0<y 

and define 5 to be 1 if f > 0, if f < 0, and any integer if f = 0. By the 
definition of particle current fl3T|) - fl33|l . we have that, for any C ^ X, 

Mr(C,c.,a;+) := ^f-^'^^'iCu) 







-v^^.{C,^,^^)ix + 5)]ds (55) 
Er\Cu;,u;^) := exp {^^f (-")+^(C, .;) 

-(e^-1) f r{T,?f,-{C,u,u+))ds 
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^0 

- £ \v\ (^e-'^-(^n-(C,^'-^)(^+') - l) ds^ (56) 
exp{9Mr{C,uJ,u;+) 

Jo 





i 



^1 |^g-sgn^■^)JPT;^_^(,,^^,aJ■ j^x-hoj _ 

+sgn(w)6'r/^„(C,u;,u;+)(x + 5)) ds] 

are martingales under P ® P'^, with respective means and 1. Notice that 
rjs- and rj^ can be replaced with each other in the above martingales, because, 
by the graphical construction of Section [2711 s ^— ?7s(x) is P(g>P^-a.s. locally 
piecewise constant for every C ^ X and x G Z. It follows from ( l56l) that 

E (e^*^*"") < e^*(-"'^'~i-^^'l^l) (57) 

for any G X, where expectation is w.r.t. P ® P"*", and the constant C 
depends only on p(.), 6(.) and f but not on ^. Cramer's inequality and (IFTI) 
imply the large deviation bound 

P 8) P+({|Mf ''^l > y}) < 2e-*^''(^) (58) 

for any C ^ X and y > 0, with the rate function 

Icfe) = flog (l + log (l + ^)] 

Because of the independence property fl53l) . by fl58|) . 

z7^'''®P®P+({|Mr(u7,„e[^i],,i^,^+))| >2/}) <2e-*^^(^) 
This and Borel Cantelli's lemma imply that 



limt"i|Zn [-£t,rf]r' V M't'\w^ue^pt]M^,uj+) = Q, (59) 
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(g)P+_a s In view of (1551) and fl59p . the proof of fl49l) is now reduced 
to that of the same set of inequahties with (IMI) replaced by 

t-^\Zn[-et,et]\-^ j [Tj{?li{Wat,0[(3t],atUJ,UJ+) 
" x£Z: \x\<Et 

-T^Vrf^{zUat,0lf3t],ati^,UJ^)il)]ds 

which is exactly, because of fl52|) . 

J Ifiv) - vvil)]m,,eie[f,,]^at^', u;+)idv) (60) 

with the empirical measure rrit^e defined in f l45p . By Proposition 12.11 and 
Lemma [231 every subsequential limit z/ as t ^ oo of rrit^e^O'^p^-^ ^^u',u^) is of 
the form 

for some measure 7 supported on 7^ fl [A,p]. Then the corresponding subse- 
quential limit as t ^ 00 of fl60|) is of the form 

j \G{r) — vr\'^{dr) 

as one verifies, using shift invariance of z/*", that 

j [fiv) ~ vrj{l)]v'^ {dvi) = G{r) — vr 

This concludes the proof. 

Step two: proof of ( fJOj) . Let Mi<y<t><'i;<t>i,A<pG7^, and 
r e [A, p] DTZ. We set <^ = r]^, and we define z/^''"''' as the coupling distribution 
of [t]^ , Tj'^ , r]'') ■ Note that, by the stochastic ordering property (ITSl) . 

^^'■''''{(^,^,0eX3: r/<c;<e} = l (61) 

The following limits are true for z/'*''^''^-a.e. (77,^,^). By the expression (136!) 
of A^'' and the equilibrium limit Lemma [3.11 

hm ^AT^^ o ^1 ^^,(,, u) = r{v - m) (62) 
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By attractiveness, 

liminf Vr'''oe[^,j„,(,,e,a;) > lim ^M^^^ o e[p,^^^,{,,,,u^) (63) 

t— +00 t t—*oo t 

Putting together fl62|) and (l63|) . 

liminf >r(t;-«i) (64) 

t— ►00 t 

Now, by fl36l) . Lemma [3.21 and Lemma [3.11 respectively for r and p, 

lim IaT-'- o ^|^,,^„,(,, e, uj) = (G(r) - «ir) - (G(p) - t;ip) (65) 
Subtracting ([Ml) to f[65l) . we get 

limsup V'^^ o „,(<,, e, .;) < (G(r) - vr) - {G{p) - v,p) (66) 

t— >oo t- 

By attractiveness, ( [611) and f [66l) . we have 

hm sup \Ut'''^ o ^|^t]^„t(r/, ^, ^) < hm sup \Ut'''^ o e[^i]_„i(^, ^, a;) 
t—>oo r t^oo r 

< {G{r)-vr)-{G{p)-v^p) (67) 
Using (!36|) . (!67j) . Lemma [3^ and Lemma [3?T] for p, we obtain 

limsup ° 

= hmsup o ^[^,],,,(a;') - 1^ ° + ° O'mM^^')) 

< hmsup ^A/;"'"^ o e[^,^^^,{u') + hmsup ^0- ° ^U],«t(^') 

t— ►oo r t— >oo r 

< ((G(r) - t;r) - {G{p) - v,p)) + (G(p) - v,p) 
= G{r) — vr 

for every r G [A,p] fl 7?.. Since ^ is no fonger involved in the above inequali- 
ties, we obtain a z/'^'^'-a.s. limit with respect to (?7, every r G [A, p] fl IZ. 
By continuity of G this holds outside a common exceptional set of z/'^'^- 
probability for all r G [A, p] fl 7?.. This proves fISU]) . □ 
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Proof of Corollary \ 3. 1[ 

(i). By Proposition [3l](cf. [7, (28)]), 

l-[G{K[v)) - vK{v)] = -K{v) 
weakly with respect to v. Thus, setting u = R\p, we have 

[G{u{v,l))-vu{v,l)]-[G{u{w,l))-wu{w,l)]= I u{x,l)dx (68) 

J V 

for all f , w G M. Let a < 6 in M. Setting 

W = T {tinxo]VNso{V, TlNxoWNsoi^, ^)) 

we have 

pNt{0[Nxo],NsoiV,^,^))iia,b]) = N^^ VNt{^,OlNxo],Nso^){x) 

[Na]<x<[Nb] 

Thus, by Proposition 13.21 and (l68l) . 



b/t pb 

}^^_pNtid[Nxo],Nso^')ii(^,b]) =t I u{x,l)dx= u{x,t)dx (69) 

a/t J a 



ij^'P (g) P-a.s., where the last equality follows from Proposition 13. 1[ Now (168!) 
implies that the r.h.s. of (128|) is a continuous function of {v,w), while the 
l.h.s. is a uniformly Lipschitz function of {v, w), since the number of particles 
per site is bounded. It follows that one can find a single exceptional set of 
P-probability outside which (!69l) holds simultaneously for all a, b, 
which proves the claim. 

(ii). Since r/'^''' has distribution z/'^''' under Pq, the statement follows from (i) 
with xq = Sq = 0. □ 
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3.4 Proofs of remaining lemmas 



Proof of Lemma \3.^ Let e > 0. We consider the probability space Vt' x (Z+)^ 
equipped with the product measure 

:= p^^P ® P ® 

where is the product measure on Z whose marginal at each site is Poisson 
with mean K{1 + e). A generic element of this space is denoted by (a;',x), 
with uj' = {r],C.,uj) and x ^ (Z+)^. We first prove Because of 

(that is, coupled configurations are ordered under u^'^), by the attractiveness 
property f|T3l) . we may define 

7.(0.') ■.= Vsi^,uj)-VsiTi7],0,^) (70) 

for s > 0. Therefore 7-particles represent the discrepancies between the sys- 
tem starting from ^ and the system starting from T{ri,C,)- We look for v such 
that there are no discrepancies to the right of vt, in which case the parti- 
cles there should be distributed hke ^-particles, according to the equilibrium 
measure u^. Let v > v. Because 7o(x) = for all x > 0, by the definition of 
current 

0^(e,e,^) = 0n^')+ E ^*(^) 

y>[vt] 

Therefore, to prove fl38l) . we want to obtain 

hm 7* ° 0[^t]U^')iy) = 0, ^"'^ ® P-a.s. (71) 

t— >oo ^ — f '- ' 

y>vt 

To this end, we follow the proof of pQ Proposition 5], with minor modifica- 
tions. We emphasize that even if the latter proof corresponds to a = (3 = 0, 
we will see that the arguments extend to (a, P) 7^ (0, 0). We label 7-particles 
with i?;'s and x-particles with Z 's as follows: we denote by = R^iuj'), 
resp. Zq = Zq{x)^ the position of the 70-particle, resp. x-pa-rticle, with label 
j (we take j < 0). The labelling is such that i?^ < -Rq^^ and Zg < Z^'^^ for 
all j < 0, where Rq, resp. Zq, is the position of the first 70-particle, resp. 
X-particle, to the left of (or at) site 0. By the definition of the number of 
X-particles between —n and will be eventually larger than nK. Let 

W{x) := inf {ra G Z+ : Z^ > /K] for every j < -n} 



25 



By Poisson large deviation bounds, the random variable W is P^-a.s. finite 
with exponentially decaying distribution. Since 70 (x) < K for every x G Z, 
we have < \ / K] for all j < 0, hence Zq > for every j < —W{x)- 
The dynamics of is defined by: if Rj-_ = x and, for some z > and 

ue[o,i], 

{(t, X, z, u), (t, X + z,-z,u)} n 7^ (72) 

then = +z. In other words, x-particles evolve as mutually independent 
(given their initial positions) random walks, that jump from y to y + z a.t 
rate 

piz) = {p{z)+p{-z))\\b\\^ 

for all ?/ G Z, 2; > 0. Then, since a jump for R^ from Rj-_ = x to Rj = x + z 
is possible only under ( 172|) . 

i?^ < Z^' ^ Vt > 0, Rl < Z{ (73) 

In view of ( I7T]) . (|751) . we estimate 

a:>i;t j<0 i<-H'(x) -VK(x)<j<0 

i<-i^(x) 

< Zi(a;,x) + W^(x) (74) 
where Zt := X]i<o ^{zi>vt} is a Poisson random variable with mean 

KZ, = Kil + E)J2KiYt>j) (75) 

for Yt a random walk starting at that jumps from y to y + z with rate p{z). 
Repeating P, (43)-(45)] gives that 

lim lE'^Zt/t = (76) 
if we choose v > J2z>o ^p{^)- Let 6 > 0. Since Zt is a Poisson variable, 
F',{\Zt-KZt\>6t) < 



(5t) 



[e;(z, - E^Z,)T 



4 



m 

< ^ ^ ^ ^ — (77) 

{Sty ^ ^ 
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Therefore, by Borel Cantelli lemma 

lim (Zt - KZt)/t = 0, P^-a.s. (78) 



Since P is invariant by 6[i3t],at, Zt{9[f3t],at^^ , x) has the same distribution as 
Zt{uj,x) under P^. Thus fl75|) - (178|l still hold with Zt{9[i3t],at^ , x) instead of 
Zt(u;,x),and 

lim t-^Zt{e[pt]^^too, x) = 0, P;-a.s. (79) 

t — >oo 

Because the random variable W in fITil) does not depend on u', fl79l) and fl7^ 
imply fITT]) . This concludes the proof of 



If we now define 7^ as 

7,(u;') := ris{T{ri, 0> ^) - Vs{v, ^) 

and replace (l)f{uj') by —(j)^{uj') (so that the current, which was rightwards, 
becomes leftwards) , then the proof of (!39l) can be obtained by repeating the 
same steps as in the previous argument. □ 



Proof of Lemma 
(i). Since 

\t]-t _ 
t J[t] 

has total variation bounded by 2/t, it is enough to prove the result for every 
subsequential limit of the sequence mn^ei^lf^n] an^'^^'^) as n ^ 00, n G N. 

Step one. We prove that every subsequential limit lies in It is enough 
to show that, for every open neighborhood O of 2^, with u^'^ ® P ® P^- 
probability one, mn,e lies in O for sufficiently large n. One can see from (H5!) 
and (HZD that 

where, for fixed t, ^* is the process defined by 
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with the configuration Wat defined in and satisfying (jHSD- Hence under 
^A,p^p^p+^ -g ^ Markov process with generator and initial distribution 
v^'P independent of t. By Lemma I3.4[ 

hmSUpn-Mogl?^'^® P ®P+ {mn,e{0[pn],an^\^^) ^ O) < 
n^oo 

Now Borel-CanteUi's lemma implies that, a.s. with respect to z/'^'''(g>P®P^, 
Step two. We prove that every subsequential limit lies in M.\^p- Since 

for v^'P- a.e. [1],^,), fl52p and the monotonicity property (1131) imply 

for all t > 0. By dSD and 



^0 

By dSD-dHl) and (00]) -dH]), 



Thus 



/at+t 
\Zn[-et,et]\-' r,.(5^ 



u^'P (g) P (g) P+- a.s. for all t > 0, where 

^"^(0;+) := Z n [[/3t] - sgn (t;)Ar,t(cu+) - et, [fit] - sgn (t;)Ar,t(cu+) + et] 



28 



We now argue that It and respectively converge a.s. to u and with 
respect to z/'^'^ ® P ® P^. Let us consider for instance It- Let Lf denote the 
measure defined as k but with Xt{uj~^) replaced by 

yt:=Zn [{p -av- e)t, {p - av + e)t] 

By the strong law of large numbers for Poisson processes, there exists a sub- 
set C C with P'^-probability one, on which sgn{v)Nt{uj^)/t ^ f as 
t ^ oo. The total variation of It — Lt is bounded by {2et + |A'(A3^j|, 
where A denotes the symmetric difference of two sets. Hence, for uj^ G C, 
It — Lt converges to in total variation. We are thus reduced to proving a.s. 
convergence of Lt to u^. Under V^'^^P (S>P^, Tfif], uj, uj'^) is a Feller process 
with generator Ly and initial distribution z/^. Thanks to (ITS!) and (H3!) . we 
can apply Proposition 12. 3[ or more precisely, its extended form (!22l) . This 
implies convergence of Lt. 

Step three. We prove that every subsequential limit lies in S. To this end we 
simply note that the measure 

Timn,e-mn,£ = \[-en, ETl] D \ ^ T^rrin 

ya;eZn(-en+l,en+l] 

- ^ r^rrin (81) 

x£l,n[—en,en) I 

has total variation bounded by 2 |[— eri, e'n\ fl Z|~^ = 0(l/n). Letting n — ^ oo 
in ( [HT]) shows that Txvn = m for any subsequential limit m of m^ g. 

(ii). Proposition 12.11 implies u = J u^'~f{dr) with 7 supported on TZ. Let 
[A', p'] C TZ denote the support of 7, and assume for instance that A' < A. 
Choose some A" G (A', A). By Proposition 12.21 the random variable 

I 

M{r]) := lim(2/ + 1)"^ V t]{x) 

x=—l 

is defined u'^-a.s. for every r G 7^, and thus also u-a.s. It is a nondecreasing 
function of rj. Thus, <v implies 

z/(M < A") < z/\M < A") = 
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where the last equahty follows from Proposition 12.21 hence a contradiction. 
Similarly p' > p implies a contradiction. Thus 7 is supported on 7?. fl [A,p]. 
□ 



Proof of Lemma \3.4 



(i) . Nonnegativity follows from taking / = in (jlSD- As a supremum of 
continuous functions, is lower semicontinuous : indeed, because the num- 
ber of particles per site is bounded, each local function is continuous and 
bounded and so is L^,(e-^)/e-^, hence the functional defined on P(X) by 

= - 1 ^d^, (82) 

is continuous. The inclusion C V~^{fS) holds because of 

Lvi}ogg) < L^g/g 

which follows from the elementary inequality log h — log a < {b — a)/a, and 
the fact that / Ly(\ogg)dp = ii p E ly. We eventually prove the reverse 
inclusion V~^{0) C X„. Fix a local test function /. li p E V~^{0), we must 
have 

I{t) := / U p > 0, yteR (83) 

As / is local and the space {0, . . . , K} is finite, integrability conditions are 
satisfied to differentiate I(t) in fl83|) under the integral. Since /(O) = 0, 
equation fl83l) implies that I{t) has a minimum at t = 0, hence 

= — — = / — — — dp = / LJdp 
dt |t=o J dt e-J |i=o J 

and thus p & 1^, since this is true for any local function. 

(ii) . Since 0/ is continuous, by [24, Appendix 2, Lemma 3.3], it is enough to 
prove that 

limsuprMogP Ut,e{i) e O) < inf sup-0/(/i) (84) 

t-+00 ^ ^ / local ^gO 
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for every open subset O C 'P(X). Let / be a local test function on X, and 

set 

oo 

/(t,r/) := |Zn [-£t,£t]r^ ^ =^l[ne-\{n+l)e-^){t)fn{r]) 

x&n[-£t,et] n=0 

where 

n 

/„(r7) := (2n + l)-i J]r./(ry) 

For each n G Z"*", 

M/'" := exp iMt) - - [ e-^~"L4e^~"]{l-)ds\ , t > ne"^ 



is a mean 1 martingale under P with respect to the cr-field Qt generated by 
(6, s<t) (cf. [211 Section 7 of Appendix 1]). It follows that M/ defined for 
t > by 

n 

Mi := H Mg-i)_ M/'", t G [ne-\ (n + 
fc=i 

(where the product is 1 for n = 0) is a mean 1 ^t-martingale under P. Thus 
we can define a probability measure P-'^ on Qt by dPf/dP = M/. A simple 
computation shows that 



M/ = exp /(t,ei)-/(0,eo)- / e-fL4ef]is,^,-)ds + R{ 



where 



[et] 



n=l 



-R/ = /n-l('C{ne-i)-) " fn{C,n 



^5) 



(86) 



Notice that, by the graphical construction of Section [2.11 s i-^ rjs^x) is for 
each X G Z a piecewise constant function whose jumps occur at (random) 
times which are a subset of some Poisson process. Thus s~ in (!85|) can be 
replaced by s, and {ne~^)~ in (1861) by ne^^. The latter implies that the 
summand in (ISBI) is bounded in modulus by 4(2n + 1)~^ sup |/|. Hence 



R 



<2[l + log(£t)]sup|/| 



(87) 
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We claim (this will be established below) that, for every probability measure 
fj, on X, the mapping f ^ f L^[e^]dn (defined on the set of local functions 
/ : X — > M) is convex. Since Ly commutes with the space shift, this implies 

^<exp|-i?f + /(0,eo)-/(t,6)+i / ^(^K.(0(rf^)} 
Thus, for any open subset O of 'P(X), we have (cf. ( !82l) ) 

< exp|-i?f + 2sup|/|-Unf 0^(/i)| (88) 
Using (IHTl) and minimizing the r.h.s. of (IHHl) over local functions /, we obtain 



Proof of claim. We prove that f ^ J Ly[ef]d^ = —(j)f{fi) is convex 
for any G V(X.). Equivalently we show that, for any local functions /, g on 
X, t —(f)ff^g[fi) is convex on M. We have 



dt^ J e*/+f J e*/+9 

The above integrand is nonnegative. Indeed, for local functions ip and ip, 

Ki^^ij) - 2^L,{^i,) + ^^L,i, = Lti^p^) - 2^Lt^ (89) 

where L"^^ := L^{^i)) - ^L^ifj. For %1)>Q,^ is a Markov generator, and 
thus the r.h.s. of fl89l) is nonnegative. □ 



4 The Cauchy problem 

In Corollary 13.11 we established an almost sure hydrodynamic limit for ini- 
tial measures corresponding to the Riemann problem with 7^-valued initial 
densities. In this section we prove that this implies Theorem 12. ![ that is the 
almost sure hydrodynamic limit for any initial sequence associated with any 
measurable initial density profile (thus we add in this Section the hypothesis 
p{.) finite range, which was not necessary for Riemann a.s. hydrodynamics). 
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This passage is inspired by Glimm's scheme, a well-known procedure in the 
theory of hyperbolic conservation laws, by which one constructs general en- 
tropy solutions using only Riemann solutions (see e.g. [371 Chapter 5]). In 
[Zt Section 5], we undertook such a derivation for convergence in probabil- 
ity. In the present context of almost sure convergence, new error analysis is 
necessary. In particular, we have to do an explicit time discretization (vs. 
the "instantaneous limit" of [HI Section 3, Theorem 3.2] or [?', Section 5] 
for the analogue of fllOOl) below), we need estimates uniform in time (Lemma 
14.21) . and each approximation step requires a control with exponential bounds 
(Proposition 14.21 and Lemma 14.31) . 

4.1 Preliminary results 

For two measures G A^^(M) with compact support, we define 



When a or /3 is of the form u{.)dx for «(.) G L°°(R) with compact support, 
we simply write u in fl9UI) instead of u{.)dx. A connection between A and 
vague convergence is given by the following technical lemma, whose proof is 
left to the reader. 

Lemma 4.1 (i) Let {a^)N be a M.^ (W) -valued sequence supported on a com- 
mon compact subset ofR, and u{.) G L°°(R). The following statements are 
equivalent: (a) u{.)dx as N oo; (h) A(a;^,M(.)) ^ Q as N oo. 

(a) Let (a;^(.))iv be a sequence of A4^(R) -valued functions : T — ^ 
A1"'"(]R), where T is an arbitrary set, such that the measures Q;^(t) are 
supported on a common compact subset of M. Assume that, for some a : 
[0, +oo) Ai^(M.), A{a^ (t) , a{t)) converges to uniformly on T. Then 
<y^{.) converges to a{.) uniformly on T. 

The following proposition is a collection of results on entropy solutions. We 
first recall two definitions. A sequence (m„, n G N) of Borel measurable 
functions on M is said to converge to u in Lj'jjj,(M) if and only if 



A(a, /3) := sup |a((— oo, x]) — /5((— oo, x])| 



(90) 
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for every bounded interval J C R. The variation of a function u{.) on an 
interval / C M is defined by 

{n-l 
\u{xi+i) - u{xi) I : e M, Xo, . . . , x„ G /, Xo < ■ ■ ■ < Xn 
i=0 

We shall simply write TV for TV^. We say that u = u{., .) defined on Rx R"*"* 
has locally bounded space variation if for every bounded space interval / C R 
and bounded time interval J C R"*"* 

sup TV/['u(., t)] < +00 

Proposition 4.1 

0) Let u{., .) be the entropy solution to ^ with Cauchy datum uq G L°°(R). 
Then the mapping t 1— >• u{.,t) lies in C°([0, +00), Ll^^{R)). 

1) If Uq{.) is a.e. IZ-valued, then so is the corresponding entropy solution 
u{.,t) to ^ at later times. 

a) If Uq{.) has finite variation, that is TVMg(.) < +00, then so does u\.,t) 
for every t > 0, and TV-u*(.,t) < TVmo(.)- 

Hi) Finite propagation property: Assume (i G {l,2}j is the entropy 

solution to l^j with Cauchy data Uq{.). Let V = ||G'||oo := supp|G"(p)|. 
Then: (a) for every x < y and <t < {y — x)/2V , 

ry-Vt r-y 

/ \u^{z,t)-u^{z,t)\dz< I \ul{z) -ul{z)\dz (91) 

Jx+Vt Jx 

In particular, ifu^ is supported (resp. coincides with Uq) in [—R, R] for some 
R> 0, u^{., t) is supported (resp. coincides with u'^{.,t) ) in [—R~Vt, R+Vt]. 
(b) U lKui{z)dz < +00, 

A{u\.,t),uX,t))<A{ul{.),ul{.)) (92) 

iv) Let xq = —00 < a;i < • • ■ < a;„ < Xn+i = +00 and e := mino<fc<n(a;fc+i — 
Xk). Denote by Uo{.) the piecewise constant profile with value on Ik '■ = 
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(xfc, Xfe+i). Then, fort < e/{2V), the entropy solution u{.,t) to ^ with 
Cauchy datum Uq{.) is given by 

U{x,t) = Rr^_^^rk{x - Xk,t), Mx G (Xfc_i + Vt,Xk+l - Vt) 

Properties o), ii) and iii) are standard, see [2ZIEZIEB]- Properties i) and iv) 
are respectively [H Lemma 5.3] and P, Lemma 3.4]. The latter states that 
the entropy solution starting from a piecewise constant profile can be con- 
structed at small times as a superposition of successive non-interacting Rie- 
mann waves. This is a consequence of iii). Note that the whole space is indeed 
covered by the definition of t) in iv), since we have Xk+i — Vt>Xk + Vt 
for t < e/{2V). 

The next lemma improves [71 Lemma 5.5] by deriving an approximation 
uniform in time. 

Lemma 4.2 Assume uo{.) is a.e. TZ-valued, has hounded support and finite 
variation, and let {x,t) u{x,t) be the entropy solution to ^ with Cauchy 
datum uq{.). For every e > 0, letVe be the set of piecewise constant TZ-valued 
functions on M with compact support and step lengths at least e, and set 

6,{t) :=£-Mnf{A(n(.),n(.,t)) : u{.) G 

Then there is a sequence | as n oo such that 6s„ converges to 
uniformly on any bounded subset o/R+. 



Proof of Lemma 4-2 We first argue that, for every e > 0, is a continuous 
function. Indeed, by Proposition 14.11 iii), a) for every T > 0, there exists a 
bounded set Kt C M such that the support of u{.,t) is contained in Kt for 
every t G [0,T]. Since 

A{v,w) < / \v{x) — w{x)\ dx (93) 
Jr 



for v,w E L°°(M) with compact support, it follows by Proposition 14. 11 o) and 
Lemma [4. H (i) that 

limA(u(.,s),M(.,t)) = (94) 
for every t > 0. This and the inequality 

Mt) - 6,{s)\ < e-'A{u{.,s),u{.,t)) 
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imply continuity of 5^. By Proposition 14. i) and ii), u{.,t) has bounded, fi- 
nite space variation, and is 7^- valued. Hence, by [H Lemma 5.5], for any given 
5 > 0, for e > small enough, there exists an approximation e Ve 

of u(.,t) with A {u^'\.,t),u{.,t)) < eS. This implies d^it) as e ^ for 
every t > 0. Let T be some countable dense subset of [0, +oo). By the diag- 
onal extraction procedure, we can find a sequence En I such that ^^^(t) i 
for each t E T. By continuity of 5^ we also have that 5e„(t) i for every 
t G [0, +oo). Dini's theorem implies that 6^^ converges uniformly to on 
every bounded subset of [0, +oo). □ 

We now quote P Proposition 3.1], which yields that A is an "almost" non- 
increasing functional for two coupled particle systems: 

Proposition 4.2 Assume p{.) is finite range. Then there exist constants 
C > and c > 0, depending only on b{., .) and p{-), such that the following 
holds. For every N G N, (r^c^o) e with |r7o| + |^o| := Y.xez['no{x)+^o{x)] < 
+00, and every 7 > 0, the event 

Vt > : A(a^(r/,(r/o,^)),a^(^*(eo,^))) < A(a^(r^o), «^(eo)) + 7 (95) 
has F -probability at least 1 — C(|?7o| + |'Co|)e~'^^'^. 

We finally recall the finite propagation property at particle level (see [TJ 
Lemma 5.2]), which is a microscopic analogue of Proposition 14. iii). 

Lemma 4.3 There exist constants v and C, depending only on b{.,.) and 
p{.), such that the following holds. For any x,y E Z, any (r^c^o) ^ X^; CLnd 
any < t < {y — x)/{2v): if rjo and ^0 coincide on the site interval [x,y], 
then with F -probability at least 1 — e"*"*, rjs{rjQ,u) and r]s{^Q,u!) coincide on 
the site interval [x -\- vt,y — vt] fl Z for every s G [0, t]. 

Remark 4.1 The time uniformity in Proposition \4.2\ and Lemma \4.3\ does 
not appear in the original statements (repectively, |3, Proposition 3.1] and 
1^ Lemma 5.2]), but follows in each case from the proof. 

4.2 Proof of Theorem [231 
4.2.1 Simplified initial conditions 

We will first prove Theorem 12.11 under the simplifying assumptions: 

Mo is a.e. 7?.-valued (96) 
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TVmo < +00 (97) 
and there exists R > (independent of A^) such that 

supp Mo C [-R, R] (98) 

VA^ G N, Po {Voi^) = whenever x G Z, |a;| > RN) = 1 (99) 

The essential part of the work (that is, the approximation scheme) is con- 
tained here, and the proof under general assumptions will follow in Subsection 
14.2.21 by approximation arguments. 

Assumption (jHSD implies by Proposition 14. 11 i), that u{.,t) is 7^-valued, (jHT]) 
by Proposition 14. ii), that u{.,t) has finite variation for t > 0, and (198!) by 
Proposition 14.11 iii), that u{.,t) is supported on [—{R + Vt),R + Vt]. In the 
sequel we abbreviate, for {ljo,lj) G Q 

vf = VtiVo K),^) 

We consider the random process on Q 

A^(t) :=A(a^(r/^,),^.,t)) 

By initial assumption ([7j) and (i) of Lemma 14. 11 A^(0) converges to 0, Pg-a.s. 
Fix an arbitrary time T > 0. We are going to prove that 

lim sup A^(t) = 0, P-a.s. (100) 

^^°°te[o,T] 

Then one can find a set of probability one on which this holds simultaneously 
for all T > 0. Theorem 12.11 follows from (ii) of Lemma [4. 1[ 

Let e = En he given by Lemma 14. 2^ and 6 = 6n = 2 supjgp j.] 5e„(t), so 
that — ^ as n ^ 00. In the sequel, for notational simplicity, we omit 
mention of n. We fix a time discretization step 

e' = 8 mm{{2v)-\ (2^)"^) (101) 

where V and v are the constants defined in Proposition 14.11 and Lemma 14.31 
Let tfc = ke', for A; < /C := [T/e'], tjc+i = T. The main step to derive (11001) 
is to obtain a time discretized version of it, namely 
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Lemma 4.4 

limsup sup [A^(4+i) - A^(tfc)] < 3fe, P-a.s. 

Af^oo k=0,...,K.-l 

The second, more technical step, will be to fill in the gaps between discretized 
times by a uniform estimate for the time modulus of continuity, that is 

Lemma 4.5 

lim limsup sup sup A la'^ (rj^^) , {r]^^ )1 =0, P-a.s. 

e=e„^0 Ar_^oo k=0,...,K. te[tk,tk+i] 

By o) of Proposition 14. t ^ u{t,-) is uniformly continuous from [0,T] to 
Lj^g^(M). This and fl93p imply an analogue of Lemma 14.51 at the level of 
entropy solutions, namely 

lim sup sup A{u{.,t),u{.,tk)) = (102) 

Then f llOOp follows from Lemma I4.4[ Lemma 14.51 and (11021) . 

Proof of Lemma \4.4\ The method is to approximate u{.,tk) by an 7^-valued 
step function Vk{-), to associate to the profile Vk{.) a sequence of configura- 
tions {C^''')n (in the sense ( 11131) below), to use Riemann hydrodynamics for 
these approximations from time tk up to time t^+i, and to show that approx- 
imated systems at time t^+i are close enough to the original ones, both at 
microscopic and macroscopic levels. 

Let Vk{-) be an approximation of u{.,tk) given by Lemma [4. 2 [ so that 

A{u{.,tk),Vki.))<5e, A; = 0,...,/C-1 (103) 

We write Vk{.) as 

Ik 

1=0 

where -oo = Xk,o < Xk,i < . . . < Xk,i^ < Xk,if^+i = +oo, rk,i G TZ, rk,o = 
rk,if, = 0, and, for 1 < / < k, 

Xk,i - Xk,i-i > e (105) 
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For tk <t < tk+i, we denote by Vk{.,t) the entropy solution to at time t 
with Cauchy datum Vk{-)- For I = 1, . . . ,1^, define on Q the configurations 

e^''KuJo,uj) and ^""^'{000,00) by 

e^''U,u.)ix) := ( ^-^-f^'V^'^^'^^f ) ' < [^"^''j (106) 

e'^''=(a;o,a;)(x) := rimM'"'i^o),^)ix), if [iVx^,;] < x < [Nxk,i+i] (107) 
so that ^^''^(co'o, t^) has finitely many particles, and 

^N,k^^^ = ^N^k>^(x), if [iVa;fc,,_i] <x< [Nxk,i+i] (108) 

Moreover by the commutation property (jH]), 

^[iVx,,]e'^''''(^o,^) (109) 

Evolutions from fll06p - fll07l) are denoted by 

ef''^'^(^o,^) = Vt{^''''\uo,u^),eo,Nt,u^) (110) 
so that by the Markov property and fll09p . we have 

nNx^,i]^^'^'\u;o,u;) (111) 

We claim that 

lim a^{r]mk{'n'~''''{^o),^)){dx) = rk,idx, P-a.s. (112) 

N^oo 



For k = 0, this follows from Proposition 12.21 since //'"''■'(tuo) ~ For 
= 1,...,/C — 1, this follows from Corollary 13.11 with X = p = r^^i and 
So = Xq = 0. Indeed, on the one hand we have i?rfc ,,rfe ,(•' •) = on the 
other hand, if u' = {r],C,,uj) ~ z/^fe.*'^'*:.' ^ we have f] = ^ a.s., so that 
(a;') = a^'ivNtdv,^)) a.s., with (r/,a;) ~ u^^-^ ®1P. 
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By fllQ7p . for every continuous function ■?/' on M with compact support, 



1=1 -^^ 



){x)a^{VmM'"M,oo)){dx) + 0{l/N) 



N^oo 



i'{x)vk{x)dx, P-a.s. 



where the convergence follows from (11121) and (11041) . Hence, 
lim a^{^^'''{uJo,uj)){dx) = Vk{.)dx, P-a.s. 



;ii3) 



that is, ^^''^ is a microscopic version of Vk{ )- For /c = 0, . . . , /C — 1, we write 
(remember that e' = t^+i — t^) 

A^(4+0-A^(4) < A[a^(ryi^,^^J,a^(e;j^)l-A^(4) 

a^(eS, .;.(., .') 



+ A 

+ A{vk{.,e'),u{.,tk+i)) 
By (I103P and iii), b) of Proposition 14. H 

A(t;fc(., tfe+i - tfc), m(., < A{vk{.),u{.,tk)) < Se 



(114) 



;il5) 



is a bound for the third term on the r.h.s. of (I114p . For the first term, we 
define the event 



E 



N,k 



A 



By assumption fl99l) and Proposition 14.21 

for some constant C independent of k. Thus, by Borel-Cantelli's lemma, 
there exists a random Ni{ujq,uj) such that P-a.s., E'^''' holds for all N > Ni 
and k = 0, . . . , IC — 1. On the other hand, 

A [a^(7^i^,J,a^(e^''=)] < A^(4) + A(n(.,4),t;,(.)) + A [t;,(.),a^(r''=)] 
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Thus, by ffTO . ffTT^ and (i) of Lemma SH 

limsup {A [a''{v^J,a''{e^')] - A^(4)} < 6e 

N^oo 

with P-probabihty 1. Therefore, P-a.s., for = 0, . . . , /C — 1, 

-A^(tfc)} <2fc (116) 

is a bound for the first term on the r.h.s. of (11141) . We finally bound the 
second term on the r.h.s. of (11141) . For = 0, . . . , /C — 1 and / = 1, . . . ,lk, 
by the respective definitions ([26]), ([29]) of 6', (3^, and ffT06|) . ffTTOj) . ffTTTj) we 
have 

(7--[jvx,,]/iv)a^ (^^^t'\uJo,u)^ = ((^[7Va.fe,,])^i^f ''(^o, ^)) 
= /5i^,o^[^,^^^j_^,^(r^'^'=.-^(u;o),r7'-'=^'(^o),^) (117) 
for any t > 0. This and Corollary 13. II imply 

\im U^f'^) = Rr^ r^ {. - Xk,i,t)dx, P-a.s. (118) 

Let us consider the events 

pN,k,i |^jv,fc(^^ ^ ej^'(a;), WxeZn [N{xk,i-i + ve'), N{xk,i+i - ve'))] 
By ffTUg]) . the definition ffTUB of e', ffTUS]) and Lemma SSI we have 

win- F^''^'') < e-^^^' 

Thus there exists a random A'"2(ix;o, cu) such that P-a.s., F^'^'^ holds for every 
N > iV2, k = 0, . . . , /C - 1 and / = I, . . . Jk- This combined with ffTTSD 
implies that P-a.s., the restriction of a^{C,^'^,) to (a;^,;-! + ve\xk^i+i ~ ve') 
converges as iV — >• oo to the restriction of Rr^ i-i,rk X- ~ Xk,i,£')dx. By (llOlj) 
and iv) of Proposition 14. H this induces 

lim a^{C,^'!^,) = Vk{.,e')dx, P-a.s. 

which, by Lemma 14.11 implies that the second term on the r.h.s. of (11141) 
converges P-a.s. to as ^ oo. Together with (11151) and (11161) . this yields 
Lemma 14.41 □ 



lim sup 



A 



N/ N \ N / ,-N,k\ 
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Proof of Lemma \4-5\ We label ?7-particles increasingly from left to right 
at each time Ntk, denoting their positions by {R'^'^)i^i, where / is a finite 
set whose cardinal |/|, of order 0{N) by assumption fl99l) . is the number of 
particles in the system. For simplicity we omit the dependence of the labelling 
on in the notation. The position of particle i at time 6 G [Ntk, Nt^^i] is 
denoted by Rq'\ Let for any s, t G [tk,tk+i], 



A. 



A(a^(r/^J,«^(r/^J) = iV 



Nf^N 



r-1 



sup 5^ [vm{y)-VNs{y)\ 



X&1 



y<x 



Let 2; G Z be a point at which the supremum above is attained. We can 
suppose without loss of generality that 

NA,, = Y,vUy)-Y.'^m{y)- 

y<z y<z 

Therefore to the left of z at time Ns there are at least NA^^t more particles 
than at time Nt. Let Ig and It be the subsets of / which label the particles to 
the left of or at z at times Ns and Nt respectively. Then we have — |/t| > 
NAs^t which implies \Is\It\ > NAs,t- Now if i G Is\It, 



> z since i ^ h 
R^/ < z since i G /, 



(119) 
(120) 



By flll9p . since we have at most K particles per site, maxjg/^^/^ R^\ > z + 
K-^NAs,t. This implies maXi^j^\i^{R^^^ - i^^'J > R-^NA^^t by thus 



Ksup(i?^*, 



i€l 



R%) > NA,, 



and we conclude that 



A(«^(r/^J,«^(r/i^,))<iriV-Sup 



Proceeding as in the proof of Lemma [3.21 it is possible to construct processes 
Q^'^ and S'''^ on the time interval [Ntk, Ntk+i] such that 



Nt 
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for t G with: S"^'* (resp. Q'^'*) is a Markov process on Z starting 

from at time Nt^, that jumps from x to x + 2; at rate p{z) \ \b\\^ only for 
z > (resp. only for z < 0). Therefore, 

Pjsup sup A{a''{v^,),a''{v^J)>Ce 
\ k te[tk,tk+i] 



k iel 

^ EE r 

+ EE r(s 

A; i£l 



sup 



> CNe 



-Q 



k,i 

N{tk+i-tk) 



k,i 



N{tk+i-tk) 



> 



> CNe 



CNe] 



;i21) 
;i22) 



Since p{.) has finite first moment, by large deviation bounds for random 
walks, the constant C can be chosen large enough such that the probabilities 
in fll2ip - fll22l) are smaller than e"*" for some constant C (recall that 
tk+i —tk=e'isa. multiple of e). By Borel Cantelli's lemma we conclude that 

limsup sup sup A (a^(?7j^J, ^^(r^j^j^)) < Ce 

N^oo k=0,...,JCte[tk,tk+i] 

for e = en on a. set of probability one, which can be chosen common to all 
(the countably many) values of n eN. On this set we thus have Lemma 14.51 
□ 



4.2.2 General case 

We will relax assumptions fl96|) - fl99|) in two steps. 

Step one: compact support only. We prove Theorem 12.11 when the addi- 
tional assumptions flUS]) - flU^ are maintained, but fl^ - flU7|) are relaxed. Let 
T > 0. By approximating the initial profile by 7^- valued ones, we define a se- 
quence ('Uo)ngN of [0, ii"]-valued functions satisfying (p8l) - (!99l) and (p6l) - (!971) 
for fixed n, such that 

lim A(m^,mo) = (123) 

n^oo 

and a family of (deterministic) particle configurations (^7o '^)neN,AreN satisfy- 
ing (1^ for fixed n, such that 

hm A(«^(r^o"'''),<)=0 (124) 
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for each n G N. Indeed, let us partition [—R, R] into finitely many intervals 
In,k of length at most 5„ — > 0, and set 



(•^n,fc i^n,k~^ Pn,k^n,k /-^) 



where ln,k denotes the length of In,k^ Xn,k its left extremity, and pn,k the 
mean value of Uq on /„ fc. Then Uq has the same mean value as Mq on fc, 
hence A(uq, uq) < K5n- Then we define a sequence of particle configurations 
associated to Mq by 



We denote by u'^{x,t) the entropy solution to IQ at time t starting from 
Cauchy datum Uq, and by r^"'^ := r^f (t^q'^, tu) the evolved particle configura- 
tion starting from ^^g'^. By triangle inequality for A, 



A(a^(ry^,),^.,t)) < A(a^«),a^(ry;^f; 



+ A(«^(r/;^f),n"(.,t)) 

+ AK(.,t),M(.,t)) (125) 



We have by iii), b) of Proposition 14.11 

AK(.,t),n(.,t))<AK,no) (126) 

Further, by Subsection 14.2.11 fll24|) implies the analogue of fllOOl) for r^*^'^, 
that is 

lim sup A fa^(r7;^f),u"(.,t)) = P-a.s. (127) 
On the other hand, 

A (a^(ry^J, a^(r^-f )) = A (a^(r^o^), a^(ryo"''^)) + (128) 



where, by Proposition 14.21 Pj^^" is a random variable which satisfies 

P (sup P^f > 7 ) < C'Ne-'^^, V7 > 
\t>o J 
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for some constant C > independent of n. Applying Borel-Cantelli's lemma 
to a vanishing sequence of values of 7, 

lim sup r^f = (129) 

N-*oo f>o 

P-a.s.. Furthermore, 

A(a^(r/o^),a^(r/o"''')) < A {a'' {v^),Uo) + A iuo,u-) 

+A(«o",a^(r/o"''')) (130) 
By (CD, (i) of Lemma O and ffTMD-ffTaOD. 

limsup sup A [a^{r]^^),u{.,t)) <2A{uq,Uo) 

N-*oo te[0,T] 

on a subset of Q with P-probability one, which can be chosen to be the same 
for all (countably many) values of n G N and T > 0. The conclusion of 
Theorem 12.11 then follows from (11231) and (ii) of Lemma 14.11 

Step two: general case. We now finally relax assumptions fl98|) - fl99l) . thanks 
to the finite propagation property (both at microscopic and macroscopic lev- 
els). Consider uq and tiq as in the statement of Theorem 12.11 without any 
restriction. Let w = max(y, f), where V and v are the constants given 
respectively in Proposition 14.11 and Lemma 14. 3[ For n G N, we set 

Uq := Mol[-n,n.], Vo'^{x) = Vo {x)lzn[- Nn,Nn]{x) 

By Lemma [4.31 and Borel-Cantelli's lemma, P-a.s. for large enough A^, 

r^^^(x) = Vn^{x), Vt < n/{2w), Vx G [-A^n/2, A^n/2] n Z 

By the previous step, for each n G N, P-a.s., {rf^^) converges to m"(., t)dx 
as N ^ 00, uniformly on bounded times intervals. By iii), (a) of Proposition 
El for every t < n/{2w), M"(.,t) = u{.,t) on [-Nn/2, Nn/2]. Thus, for 
every continuous functions : M ^ M supported on [— A^n/2, A^r2/2], there 
is an event of P-probability one on which 

Jr 
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uniformly on the time interval [0, n/(2w)]. This event can be chosen to be the 
same for all values of n and for a countable set of continuous functions with 
compact support that is convergence determining for the vague topology. 
This establishes the result. 



A Proof of Corollary 12.1 



Let /if^ denote the distribution at time t of a Markov process with generator 
([T]). Assume {vi){dx) converges in /i^-probability to UQ{.)dx, that is, for 
all e > and every continuous function on M with compact support, 



N 



lim yUp 



N 



T] : 



ilj{x)a {ri){dx) — / ilj{x)uo{x)dx 



> e 







Then for every t > 0, a {ri){dx) converges in /xj^j-probability to u{.,t)dx. 
This weak law follows from the strong law in Theorem 12.11 Indeed, by Sko- 
rokhod's representation theorem, we can find a probability space [Qq, JFq, Pq) 
and a sequence [7]^)^ of X- valued random variables on Qq such that rj^ has 
distribution /i^, and {rjQ){dx) converges Po-a.s. to UQ{.)dx. 



B Remarks on subadditivity 

As outlined below, it would be possible to establish fl37j) in the particular case 
/3 = a = by using the subadditive ergodic theorem as in [I, Proposition 3]. 
However we cannot use this approach when a) ^ (0, 0). 



Let us introduce 



Xo,„(^') := ct>li,{J) - ^li,{r^,uj) (131) 



■ , 0,n-m["rn,m/v^ , 

then Xm,n is the same as defined in equation [H (27)] and, by [H p. 226], it 
satisfies the superadditivity property 

Xo^n > Xo^rn + Xm,n (132) 

(superadditivity is obtained here rather than subadditivity in [T], because we 
have X < p instead of A > p). We point out that the proof of fll32p in [T] uses 
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only attractiveness and the fact that we start with rj < ^, but not the choice 
of the distribution of {ri,C,)- It can thus be generahzed from the asymmetric 
exclusion process to our setting. Let us now assume that the probability 
measure on Q' is u^''' ® P. We can proceed as in [Ij. Indeed, because u^'^ is 
invariant for the coupled process, fl26|) implies that z/'^''' ® P is invariant by 
the shift 6'^ ^. By (^^, fll32p is true jj^'P ®]P-a.s. This and Poisson bounds on 
the expectation of Xo,„ imply, by Kingman's subadditive ergodic theorem, 
that n~^Xo^n{^') converges u^'^ ® P-a.s. On the other hand, n~^(f'^^^{ri,uj) 
converges u^'^ ® P-a.s. by (13.11) below. Hence, 

u^'f" ® P a.s., 3 lim n-V^^(^') (133) 

The limit in (I133P can then be identified using the hydrodynamic limit of [7] , 
in the same way as [3] is used in [I]. We thus obtain a particular case of (1371) 
for /? = a = 0. However, the case (/5, a) ^ (0, 0) would require 

z/^'^ ® P a.s., 3 lim n-V^/.(^[/3n],an^') (134) 

for every /? G M and a 7^ 0. The a.s. limit (I133P only implies a limit in 
probability for the shifted current in (11340 . as the distribution of a single 
current is unchanged by the shift. In contrast the joint distribution of the 
sequence of currents may change from (11331) to (I134p : Thus we cannot simply 
derive (11341) from (I133p . On the other hand, the shifted currents Yo,n := 
X^^n ° ^'\iin\ an lougcr cujoy a super-additivity property like (11321) . so we 
cannot use the subadditive ergodic theorem to obtain (11340 . Our approach 
to obtain (I134p overcomes this difficulty by avoiding the use of subadditivity. 



C Proof of Proposition 12.3 



The main ingredient is a two-dimensional extension of Birkhoff's ergodic 
theorem: 

Proposition C.l Let {X,J-',P) be a probability space and T,S : X ^ X 
two measurable mappings such that P o = P o = P. Then, for every 
bounded J-' -measurable f : X X , the limit 



n n 

U{x) := hm - V - V fiS'T^x) (135) 
n^oo n ^ — ^ n ^ — ^ 
j=i i=i 

exists for almost every x E X with respect to P. 
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This follows from more general results established for instance in [13] or 
Chapter 6]. However we include a simpler proof adapted to our case. 



Proof of proposition ICAl By ergodic theorem there exists bounded functions 
/* and such that 



1 " 

-Y^fiS'x) := f:{x)^f.{x),P-a.s. 

i=l 

-f^/.(T^x) := /-(a;)-./„(x), P-a.5. 



(136) 
(137) 



By f ll36p and Egorov's theorem, there is a sequence of subsets G J-', k E N, 
such that z/(Afc) ^ as A; — oo, and /" f^. uniformly on X\Ak. Let 



n _ n 



Then, 

\FZ{x)-fUx)\ < 
+ 



n — ' n 

j=i i=\ 



5^[/:(T^x)-/,(r^-a:)] \x\A,{T^x 
1 " 



< sup I/: - M + 2Mg\{x) + |/:,(x) - /„(x)| =: i?„,fc(x) 



where M := sup^ \ f\ and, by ergodic theorem, 

1 " • n~,oo 

9\i^) ■■= - J^'^A^CF^x) ™ gA^{x), P-a.s. 
i=i 

for some bounded, nonnegative, jF-measurable gA,, such that 



(138) 



(139) 
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By uniform convergence of on X\Ak, a.s. convergence fll37p . and fll38p . 
limsup„^^ -B„^fc(a;) < 2MgA^{x) holds P-a.s. By ( 11391) . gA^, goes to in 
L^{P) as ^ oo. Thus it has a subsequence converging to P-a.s. Letting 
/c — > oo along this subsequence concludes the proof. □ 

Proof of proposition \2.3[ 

Existence of the limit. Define the random variables Xij := /qIij^ fi'^^~^Vs)ds, 
where z,j G N and (?7s)s>o is the stationary Markov process with generator 
L and initial distribution /i. Take X = M^^^, JF the product Borel a-field, P 
the distribution of the A'-valued random variable {Xij)ij^^, T [{xij)ij^^] = 
{xij+i)ijeN, S [{xij)ij(zn] = {xi+ij)ij(zn- We have P o T"^ = P because 
(?7s)s>o is stationary, and P o = S because fi and L are invariant by r. 
Then the existence of the limit follows from Proposition IC.ll 

Identification of the limit. Let now X be the Skorokhod space of X-valued 
paths, and P = P^ the law of the Markov process with generator L and 
initial distribution /i. We consider on X the space shifts {tx)x& and time 
shifts (Tt)i>o defined as follows: if rj, = (?7s)s>o e then t^t], := {Txr]s)s>o, 
where on the r.h.s. is the spatial shift on particle configurations defined 
in Section [21 and Ttrj, := {rit+s)s>Q- What follows is a generalization of a 
standard result for one-parameter Markov processes (see e.g. [TUl Chapter 
7]). By the above existence step, we can define 

/„(r/,) := hm F^iv.) (140) 

n^oo 

P^-a.s., where 

-| pan -| " 

f:m) = — / -J2^'fivt)dt 

an Jo 

As a limit of measurable functions, is measurable. For every t > 0, 
TtF^^ — PJ^ and rP"^ — F^^ consist of space-time sums over boundary do- 
mains of order 0{n) = o(n^), hence in the limit n oo, is invariant 
by (Tt)t>Q and r := ri. To show that this implies is a P-a.s. constant 
function, we will prove that any measurable subset P of A" which is invariant 
by {Tt)t>o and r has P^j-probability or 1. Taking expectations in fll40p . the 
constant value of must he J f dfi, and Proposition 12.31 is thus established. 
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Let F G X he measurable, and invariant by {Tt)t>o and r. Set 

g{r]) = P,{F\r]o = v) =■ Pr,{F) (141) 

which is defined for /i-a.e. rj. Here, denotes the law of the Markov process 
starting from deterministic state t] We are going to prove that 

^ = or ^ = 1 (142) 

/i-a.s., which will imply P{F) = J^g{ri)fj,{dri) G {0, 1}. 

We have 

girr^) = P.,(F) = P.ir-'F) = P,{F) = g{r^) 

where the second equality follows from translation invariance fl2U]) of L (which 
implies Prri = tPti), and the third from r-invariance of F. Therefore g is /j,- 
a.s. invariant by the spatial shift r. We claim that g = Iq /x-a.s. for some 
G C X. Indeed, let J^t denote the cr-field of X generated by the mappings 
f], ^— T]s for s <t. With P^j-probability one, 

giVt) = P,XF) = P,{Ti'F\T,) = P,{F\Tt) 

where the second equality follows from Markov property, and the third from 
the Tt invariance of F. By the martingale convergence theorem, we have the 
P^-a.s. limit 

lim g{r]t) = 1f{v.) (143) 

t— >cx> 

Since rjt ^ fi for alH > 0, for every e > 0, 

P^iv. ■■ e < giVt) <l-e)=fi{v: e< g{r,) <l-e) (144) 

we conclude from fll43p - fll44l) that the law of g{rit) is Bernoulli, hence g = Iq 
/i-a.s. for some G C X. The desired conclusion f ll42p is thus equivalent to 
/i(G) G {0, 1}, which we now establish. 

First we claim that, with P^-probability one, we have g{rit) = g{vo) for all 
t > 0. Indeed, by Tj invariance of F, Markov property and definition of G, 

P,{{Vt^G}nF) = P{{^t^G}nTi'F) 

P,,{F)P,{d^) 

9{Vt)PMv.) 

{Vt(G}cX 

(145) 
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Similarly we have 

P,i{VteG}niX\F)) = (146) 

It follows from fll45p - fll46l) that {77^ G G} = F up to a set of P^-probability 
0. In other words, for every t > 0, we have g{rit) = gijio) = 1f{v.) with 
P^-probability one. 

Now, assume /i(G') ^ {0, 1}, and define the conditioned measures fiaidrj) : = 
fi{dri\ri G G) and fix\G{dri) := li{drj\rj G X\G), so that yU = ^{G)^g + 
li{X\G)iix\G- The invariance of L and with respect to space shift imply 
the same for /i^ and fix\G- A simple computation, using invariance of for 
L, and the fact that g{r]t) is a.s. constant, shows that fic and fix\G are 
also invariant for L. This contradicts the fact that G (Xfl S)e- □ 
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